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摘 要 选用 非 滑动 式 Lagrange 插值 法 和 滑动 式 Lagrange 插 值 法 对 北斗 卫星 导航 系统 Beidou 
Satellite Navigation Syetem, BDS ) 精密 星 历 进行 插值 ， 通 过 大 量 实验 对 BDS 卫星 的 插值 精 
度 与 两 种 方法 及 插值 阶 次 的 关系 进行 研究 。 实 验 结果 表明 :( 1 ) 分 别 采用 非 滑动 式 和 滑动 式 
Lagrange 插 值 时 MEO 卫 星 和 1GSO 卫 星 的 插值 误差 随 阶 次 增加 近似 呈 “U”* 形 和 “L”* 形 分 布 ， 
GEO 卫星 的 插值 误差 随 阶 次 增加 表现 出 逐渐 增 大 和 近似 平稳 的 形式 ; ( 2 ) 当 插 值 阶 次 较 低 
时 ,两 种 方法 的 插值 误差 均 会 出 现 一 定 的 规律 性 ,但 GEO 卫星 插值 误差 的 规律 性 弱 于 MEO 
卫星 和 IGSO 卫星 ;( 3 ) 各 颗 卫星 采用 两 种 方法 的 最 优 插值 阶 次 相差 均 在 1 以 内 ， 且 最 佳 插 
值 阶 次 下 ， 滑 动 式 Lagrange 插值 精度 较 非 滑动 式 可 提高 11.96%~44.01% , GEO 卫星 插值 精 
度 优 于 MEO 卫星 和 IGSO 卫星 。 
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Abstract: The non-sliding Lagrange interpolation method and sliding Lagrange interpolation 
method are used to interpolate the precise ephemeris of Beidou satellite navigation system (BDS). 
The relationship between interpolation accuracy and interpolation order of BDS satellite and two 
methods is studied. The results show that when non-sliding and sliding Lagrange interpolation are 
used respectively the interpolation errors of MEO and IGSO satellites are approximately U-shaped 
and L-shaped with the increase of the order, while the interpolation errors of GEO satellites 


increase gradually and are approximately stable with the increase of the order. When the 


interpolation order is low, the interpolation errors of the two methods have certain regularity, but 
the regularity of GEO Satellite interpolation error is weaker than that of MEO satellite and IGSO 
satellite. The difference of the optimal interpolation order of the two methods for each satellite is 
within 1, and under the best interpolation order the accuracy of sliding Lagrange interpolation is 
11.96% ~ 44.01% higher than that of non-sliding interpolation, and the interpolation accuracy of 
GEO satellite is better than that of MEO satellite and IGSO satellite. 
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0 引 


2020 年 6 月 23 日， 北斗 卫星 导航 系统 (Beidou Satellite Navigation Syetem, BDS) 全 球 
组 网 完成 , 标志 着 我 国 的 北斗 卫星 导航 系统 已 经 正式 走向 世界 并 为 全 球 服务 。 定 位 功能 是 
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区 


历 插值 精度 进 韭 行 研究 ， 对 BDS 精密 星 历 插值 的 研究 相对 较 少 。 与 GPS 卫星 星座 不 同 ， BDS 
包含 地 球 同步 融 下 轨道 卫星 ( Ce Earth Orbit，GEO)、 倾 斜 地 球 同步 轨道 卫星 
(Inclined Geosynchronous Satellite Orbit, IGSO ) 和 中 圆 地 球 轨道 卫星 (Medium Earth Orbit， 

MEO) 三 种 轨道 类 型 ， 各 类 卫星 受 插值 方法 和 插值 阶 次 的 影响 也 有 所 区 别 ， 截 止 到 目前 ， 
各 轨道 卫星 分 布 情况 见 表 1。 当 前 文献 中 ， 多 是 研究 插值 方法 的 精度 ， 并 未 对 三 类 卫星 的 插 
值 精度 情况 进行 分 析 ， 而 且 多 选用 1 颗 或 3 颗 (每 类 卫星 1 颗 ) 卫星 对 算法 进行 实验 ， 数据 
结论 的 支持 相对 单薄 "1。 

本 文 在 此 前 研究 的 基础 上 ， 选 用 滑动 式 Lagrange 插值 法 和 非 滑 动 式 Lagrange 插值 法 对 
BDS 三 类 卫星 插值 精度 进行 研究 ， 首 先 对 三 类 卫星 插值 精度 与 插值 方法 及 插值 阶 次 的 关系 
进行 研究 , 分 析 三 类 卫星 的 插值 精度 以 及 随 插值 阶 次 变化 的 情况 ; 其 次 ， 对 低 阶 次 插值 情况 
重 值 误差 的 统计 规律 进行 分 析 ; 最 后 ， 对 各 颗 实验 卫星 在 两 种 方法 下 的 最 佳 插值 阶 次 进 
行 统计 ， 并 对 最 佳 插值 阶 次 下 ， 两 种 方法 的 插值 精度 及 三 类 卫星 的 最 佳 插值 结果 进行 对 比 。 
由 于 篇 幅 限制 及 考虑 部 分 卫 星 的 数据 质量 ， 文 中 选取 PRN 号 在 C01~C37 之 间 的 32 颗 卫星 
数据 进行 实验 ， 加 粗 编号 为 选择 的 实验 卫星 。 
表 1BDS 三 类 轨道 卫星 分 布 情况 〈 加 粗 :本 文 实验 卫星 ) 


Tab.1 Distribution of three kinds of orbits of BDS satellites (bold: experimental satellites in this 
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paper) 
轨道 类 型 PRN 
GEO C01、 C02、C03、C04、C05、C18、C59、C60、C61 
IGSO C06、 C07、C08、C09、C10、C13、C16、C31、C38、C39、C40、C56 


Cll、 C12, Cl14、 C19、 C20、 C21、 C22、 C23、C24、C25、C26、C27、C28、C29、C30、C32、C33、C34、 
C35、C36、C37、C41、C42、C43、C44、C45、C46、C57、C58 
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1.1 Lagrange 插值 法 
Lagrange 插值 法 可 根据 已 知 数据 插值 出 节点 内 任意 位 置 的 数据 ， 若 有 一 系列 已 知 节 点 


和 2 及 其 对 应 的 数值 y ,yi,.….,y, ， 则 Lagrange 插值 多 项 式 可 表示 为 1: 


DT DD ee 
k=0 0 


k— Xi (Xe— Xo)... (XE— XE (XE — XE41)...(Xk— Xn) 


卫星 的 位 置信 息 中 包含 X、Y、Z 三 个 坐标 分 量 ， 则 其 在 1 时 刻 三 个 坐标 分 量 的 插值 多 
项 式 可 表示 为 : 


(一 10)...(1 一 灰 -DG 一 大 .1 一 加 ) 
X(t)= > Xx 
. 之 (不 一 10)...( 丰 一 让-D( 丰 一 友 +D (不 一 2 


(一 10)...( 一 大 -DG 一 友 +D...( 一 力 ) 
人 
的 2 (tr—f0)... (tr—te (te—terl)...(te— Xn) 


e (t—t0)...(f—tr -i)(t oti)...(t—tn) 
和 Z(1)= 》 OO 
es W 2 (tr—f0)... (te—te (te—terl)...(fe— Xn) 


(2) 


式 中 :(X(D,Y(D, 2Z0D)) 为 1 时 刻 卫 星 插值 位 置 ， tr(k = 0,1,...,n) 为 样本 点 历 元 , (X()， 


ZD，Z 有 ) 为 碌 时 刻 卫 星 位 置 。 
1.2 滑动 式 Lagrange 插值 策略 
滑动 式 Lagrange 插值 法 是 在 传统 插值 方法 的 基础 上 将 插值 区 间 视 为 一 个 可 滑动 的 “ 窗 
口 ?， 在 插值 过 程 中 令 插 值 区 间 随 待 插值 点 不 断 移 动 ， 使 待 插值 点 时 刻 位 于 插值 区 间 中 心 位 
置 呈 "29。 传 统 揪 值 方法 在 应 用 中 易 出 现 龙 格 现象 , 即 插值 区 间 端 点 附近 插值 结果 容易 出 现 震 
荡 而 导致 插值 精度 较 低 , 滑动 式 算法 仅 对 插值 区 间 中 心 位 置 进行 插值 , 能 够 避免 龙 格 现象 对 
插值 精度 的 影响 ， 进 而 提高 插值 精度 ， 但 其 计算 量 较 非 滑动 式 算 法 大 ， 在 实际 应 用 过 程 中 ， 
可 综合 考虑 插值 精度 和 效率 ， 将 滑动 式 与 非 滑 动 式 两 种 方式 组 合 使 用 。 


-> 


本 文采 用 2020 年 5 月 30 日 IGS 发 布 的 历 元 间隔 为 Smin 的 精密 星 历 作为 原始 数据 ， 并 
对 原始 数据 进行 分 析 。 依 据 表 1， 提 取 32 颗 卫 星 的 精密 星 历 ， 其 中 每 颗 卫 星 包 含 288 个 三 
维 坐标 位 置 ， 为 了 对 两 种 方法 的 插值 精度 及 插值 阶 次 进行 研究 ， 提 取 该 天 内 以 10 min 为 时 
间 间 隔 的 坐标 作为 已 知 数据 ， 对 其 余 已 知 时 刻 的 坐标 进行 插值 ， 并 与 精密 星 历 进行 对 比 , 采 
用 均 方 根 误差 〈RMS ) 作为 精度 评价 指标 。 鉴 于 实验 卫星 数目 较 多 ， 且 三 维 坐标 的 每 个 方 
向 插值 均 有 误差 ， 同 时 由 于 篇 幅 限制 ， 后 面 所 提 到 的 RMS 均 表示 卫星 位 置 综合 插值 误差 ， 
不 再 对 各 单独 方向 的 RMS 进行 统计 分 析 。 
2.1 三 类 卫星 插值 精度 与 插值 方法 及 插值 阶 次 关系 研究 

为 研究 非 滑动 式 及 滑动 式 Lagrange 插值 法 中 插值 阶 次 对 插值 精度 的 影响 情况 ， 分 别 统 
计 了 阶 次 在 4~20 阶 时 的 插值 精度 ， 结 果 见 图 1 和 图 2。 为 了 体现 三 类 卫星 在 不 同 插值 情况 
下 的 精度 量 级 ， 三 类 卫星 各 选取 两 颗 卫 星 的 插值 精度 统计 结果 列 于 表 2 和 表 3。 
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图 1 不 同 插值 阶 次 下 非 滑动 式 Lagrange 插值 法 精度 统计 


Fig.1 Accuracy statistics of non-sliding Lagrange interpolation with different interpolation orders 


RMS/mm 


图 2 不 同 插值 阶 次 下 滑动 式 Lagrange 插值 法 精度 统计 
Fig.2 Accuracy statistics of sliding Lagrange interpolation with different interpolation orders 
表 2 不 同 插值 阶 次 下 非 滑动 式 Lagrange 插值 法 精度 统计 (单位: mm) 
Tab.2 Accuracy statistics of non-sliding Lagrange interpolation with different interpolation orders 


(unit:mm) 


4 3.62 3.69 545.75 2433.90 600.65 2286.33 
9 3.64 3.37 31.17 194.27 34.47 181.68 
6 4.17 3.96 4.88 17.44 5.03 16.65 
7 4.64 4.49 5.02 5.66 4.92 4.89 
8 5.26 5.64 6.37 6.85 6.96 6.55 
9 7.44 7.20 8.62 9.52 8.41 9.17 
10 10.33 10.77 11.59 全 ,35 13.95 13.69 
11 17.56 16.24 19.26 20.38 20.04 21.71 
12 2613 21.73 26.03 34.19 30.66 30.57 
13 46.97 44.82 43.79 52.90 55.65 57.64 


14 T7733 71.10 88.98 95.38 88.43 93.39 


110.31 
245.98 
372.39 
650.75 
994.99 


2222.60 


112.19 
219.40 
358.34 
410.16 
956.14 


2204.32 


131.66 
200.47 
386.48 
673.16 
1174.84 


2199.37 


154.89 
260.19 
436.19 
752.16 
1488.41 


2747.37 


157.16 154.92 
241.98 279.74 
437.17 319.00 
931.41 688.82 
1349.69 1409.19 
2241.59 2342.38 


表 3 不 同 插值 阶 次 下 滑动 式 Lagrange 插值 法 精度 统计 单位: mm) 


Tab.3 Accuracy statistics of sliding Lagrange interpolation with different interpolation orders 


(unit:mm) 
PRN 
CO01 C05 C06 Cl11 C13 C37 
order 

4 3.02 3.12 303.99 1356.25 334.67 1274.26 
5 2.97 2.97 11.31 68.50 12.44 64.45 
6 3.00 3.03 3.49 5.32 3.15 5.16 
7 3.01 3.01 3.42 3.26 3.05 3.21 
8 3.03 3.05 3.48 3.35 3.10 3.25 
9 3.03 3.03 3.44 3.31 3.07 3.23 
10 3.04 3.06 3.49 3.36 3.12 3.27 
11 3.05 3.05 3.46 3.33 3.09 3.25 
12 3.06 3.07 3.50 3.36 3.13 3.28 
13 3.06 3.06 3.47 3.34 3.11 3.27 
14 3.07 3.08 3.51 3.37 3.14 3.29 
15 3.07 3.07 3.48 3.35 3.12 3.28 
16 3.08 3.09 3.51 3.38 3.15 3.29 
17 3.08 3.08 3.49 3.36 3.14 3.29 
18 3.09 3.10 3.52 3.38 3.16 3.30 
19 3.09 3.09 3.50 3.36 3.15 3.30 
20 3.09 3.10 3.52 3.39 3.16 3.31 
通过 分 析 可 知 : 

(1) 对 32 颗 卫 星 在 各 阶 次 下 采用 两 种 方法 插值 的 结果 对 比 , 发 现 采 用 滑动 式 Lagrange 


插值 法 进行 插值 所 获 捐 


i 值 结果 精度 均 高 于 同 阶 次 


滑动 式 Lagrange 所 


值 法 。 当 阶 次 为 4 时 ， 


两 种 方法 插值 误差 统计 的 RMS 均 比 较 大 ， 对 所 有 卫星 插值 结果 进行 统计 可 知 ， 非 滑动 式 和 
滑动 式 Lagrange 插值 结果 统计 RMS 的 平均 值 分 别 为 1629.47mm 及 908.13mm, 滑动 式 精度 


较 非 滑动 式 精度 提高 79.43%。 对 其 余 阶 次 进 


尤其 是 MEO ] 


(2)GEO | 


QD 星 和 1IGSO ] 


ED 星 ， 提 高 幅度 大 都 如 
ED 星 受 插值 方法 和 插值 阶 次 
卫星 ， 采 用 非 滑动 式 Lagrange 插值 时 ， 当 插值 阶 次 较 低 时 


9 影响 


E 50% 以 上 。 
远 小 于 MEO 


行 统计 , 滑动 式 较 非 滑动 式 插值 精度 均 有 提高 ， 


卫星 和 1IGSO 卫星 ,对 于 GEO 
，RMS 基本 都 在 4mm 以 内 ， 随 


着 阶 次 增加 ， 误 差 呈 逐 渐 增 大 的 趋势 变化 ， 当 阶 次 为 20 时 ，GEO、MEO 和 IGSO 三 类 卫星 


的 RMS 均值 分 别 为 1972.66mm、2050.71mm 和 1935.66mm，MEO ] 
计 误 差 随 阶 次 增加 均 呈 “U?” 字 形 ; 采 
增加 几乎 不 变 ，MEO 卫星 和 IGSO ] 
度 较 小 。 另 外 ， 对 GEO 卫星 的 最 佳 插 


= 


果 仅 提高 4% 左 


右 “〈 约 0.1mm) 。 


> 


滑动 式 Lagrange 插值 
在 4~5 阶 时 RMS 较 


D 星 和 IGSO 卫星 的 统 
秆 ，GEO 卫星 的 RMS 随 阶 次 的 
大 ， 随 阶 次 增加 ，RMS 变化 幅 
直 阶 次 统计 发 现 ， 最 佳 插值 阶 次 的 精度 较 4 阶 插值 结 
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(3) 非 滑 动 式 Lagrange 插值 法 受 插值 阶 次 的 影响 明显 大 于 滑动 式 Lagrange 插值 法 。 整 
体 上 看 ， 使 用 非 滑动 式 Lagrange 插值 时 ， 插 值 误差 随 阶 次 的 增加 近似 呈现 出 “U" 字 形 ， 即 插 
值 阶 次 较 低 或 者 较 高 时 ， 误 差 均 比较 大 ， 插 值 阶 次 在 7~9 附近 时 ， 误 差 较 小 ， 使 用 滑动 式 
Lagrange 插值 时 ， 插 值 误差 在 阶 次 为 4 和 5 的 时 候 ， 误 差 比 较 大 ， 随 着 阶 次 的 增加 ， 误 差 变 
化 基本 趋 于 稳定 。 
为 研究 插值 误差 的 变化 规律 ， 对 各 个 时 刻 的 插值 误差 进行 分 析 ， 顾 及 篇 幅 限 制 ， 三 类 卫 
星 各 选 一 颗 进 行 分 析 (C01、C06、C11)〉， 图 3~ 图 5 分 别 给 出 了 三 颗 卫星 不 同情 况 下 各 个 
插值 时 刻 的 误差 分 布 情况 ， 图 中 乙 轴 表 示 在 某 一 插值 历 元 下 ， 卫 星 位 置 的 综合 插值 误差 。 


x 


error/mm 
error/mm 
un 


一 20 
的 
和 一 120 100 > 
140 一 epoch 
Ca) 非 滑动 式 (b) 滑动 式 
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Fig.3 Interpolation error distribution of PRN CO1 under different interpolation orders 
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(a) 非 滑动 式 (b) 滑动 式 
图 4 不 同 插值 阶 次 下 C06 卫星 插值 误差 分 布 情况 


Fig.4 Interpolation error distribution of PRN C06 under different interpolation orders 
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(a) 非 滑动 式 (b) 滑动 式 


图 5 不 同 插值 阶 次 下 C11 卫星 插值 误差 分 布 情况 
Fig.5 Interpolation error distribution of PRN C11 under different interpolation orders 

由 图 可 知 : 采用 非 滑动 式 Lagrange 插值 时 ， 当 插值 阶 次 增加 到 10 阶 以 上 时 ,插值 过 程 
中 会 存在 部 分 插值 时 刻 误差 异常 的 情况 ， 随 着 阶 次 增加 ， 异 常 值 出 现 的 频率 相对 变 少 , 但 数 
值 却 逐渐 变 大 。 对 异常 值 进 行 分 析 ， 发 现 大 多 出 现在 插值 区 间 两 端 ， 即 异常 值 的 出 现 主 要 是 
受到 了 “ 龙 格 ” 现 象 的 影响 。 滑 动 式 Lagrange 插值 时 可 保证 插值 时 刻 位 于 区 间 中 心 ， 能 
效 避 免 “ 龙 格 ”现象 ， 故 其 出 现 异 常 值 的 概率 较 小 ， 如 图 所 示 ， 随 着 插值 阶 次 的 增加 ， 滑 动 
式 Lagrange 插值 结果 未 发 现 异 常 值 ， 误 差 变 化 均 比较 平缓。 


通过 分 析 各 阶 次 插值 误差 变化 规律 发 现 , 插值 误差 随 着 插值 历 元 的 增加 呈现 出 一 定 的 周 
期 性 ， 并 且 当 插值 阶 次 较 小 时 ， 规 律 性 表现 较为 明显 ， 对 各 颗 卫 星 三 个 方向 的 误差 作 功 率 谱 
分 析 可 知 ， 低 阶 插值 误差 的 周期 性 与 其 精密 星 历 的 周期 性 一 致 。 这 是 因为 插值 阶 次 较 小 时 ， 
重 值 所 使 用 的 原始 数据 较 少 , 各 时 刻 的 插值 结果 受 原始 数据 规律 的 影响 会 比较 大 , 随 着 插值 
阶 次 的 增加 ， 各 时 刻 的 插值 结果 受到 较 多 数据 的 作用 ， 规律 性 反而 不 会 表现 的 那么 强烈 ,这 
也 解释 了 图 中 为 何 随 着 插值 阶 次 的 增加 ， 周 期 性 反而 未 表现 出 来 。 另 外 ，GEO 卫星 插值 误 
差 的 周期 性 并 不 明显 ， 这 也 说 明了 GEO 卫星 的 插值 精度 优 于 MEO 卫星 和 IGSO 卫星 。 


2.2 两 种 方法 最 佳 插值 阶 次 精度 对 比 


同 阶 次 插值 情况 下 ， 非 滑动 式 Lagrange 插值 法 效率 较 滑 动 式 Lagrange 插值 法 高 ， 但 插 
值 精度 相对 较 低 , 为 研究 对 比 两 种 方法 的 最 佳 插值 精度 , 提取 各 颗 实 验 卫星 的 最 佳 插值 阶 次 ， 
统计 结果 见 表 4。 在 两 种 方法 各 自 的 最 佳 插值 阶 次 下 ， 对 各 颗 卫 星 的 插值 结果 进行 统计 ， 并 
以 RMS 作为 精度 评价 指标 ， 统 计 结 果 见 图 6。 
表 4 非 滑动 式 和 滑动 式 Lagrange 插值 法 最 佳 插值 阶 次 


Tab.4 The best interpolation order of the non-sliding and sliding Lagrange interpolation methods 
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PRN CO01 C02 C03 C04 C05 C06 C07 C08 C09 C10 Cll 
Non-sliding 4 5 4 4 5 6 6 6 6 6 7 
sliding 5 5 3 5 5 6 6 6 6 6 6 
PRN C12 C13 C14 C19 C20 C21 C22 C23 C24 C25 C26 
Non-sliding 7 6 7 7 7 7 7 7 7 7 7 
sliding 8 6 6 6 6 6 6 6 6 6 6 
PRN C27 C28 C29 C30 C32 C33 C34 C35 C36 C37 
Non-sliding 7 7 7 了 了 7 7 7 7 7 
sliding 6 6 6 6 6 6 7 6 6 6 
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图 6 非 滑动 式 和 滑动 式 Lagrange 插值 法 最 佳 阶 次 插值 精度 统计 
Fig.6 Statistics of the best order interpolation accuracy of the non-sliding and sliding Lagrange interpolation methods 
表 4 可知 ,各 颗 卫 星 采 用 非 滑 动 式 和 滑动 式 Lagrange 插值 法 的 最 住 插值 阶 次 均 在 4~8 
阶 之 间 ， 两 种 方法 的 最 佳 插值 阶 次 相差 均 在 1 以 内 ， 其 中 ，12.5% 的 卫星 非 滑动 式 算法 插值 
阶 次 较 低 ，59.4% 的 卫星 滑动 式 算法 插值 阶 次 较 低 ，28.1% 的 插值 阶 次 相等 。 三 类 卫星 中 ， 
GEO 卫星 的 最 佳 插值 阶 次 主要 为 4 阶 和 5 阶 ,， 低 于 IGSO 卫星 的 6 阶 ,， MEO 卫星 的 最 佳 插 
值 阶 次 相对 最 高 ， 主 要 为 6~7 阶 ， 个 别 情况 下 为 8 阶 。 
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由 图 10 可 知 ， 各 颗 卫 星 采用 最 住 插值 阶 次 插值 时 ， 滑 动 式 Lagrange 插值 精度 还 是 均 高 
于 非 滑 动 式 Lagrange 插值 ,滑动 式 Lagrange 插值 精度 较 非 滑动 式 可 提高 11.96% 至 44.01%。 
另外 ，GEO 卫星 在 插值 阶 次 更 小 的 情况 下 ， 插 值 精度 同样 优 于 MEO 卫星 和 IGSO 卫星 , 也 
说 明 GEO 卫星 的 轨道 数据 更 为 稳定 ，MEO 卫星 和 IGSO 卫星 的 插值 精度 相似 。 


3 结语 


本 文采 用 滑动 式 Lagrange 插值 法 和 非 滑 动 式 Lagrange 插值 法 对 北斗 卫星 导航 系统 的 三 
类 卫星 的 精密 星 历 进行 插值 , 对 两 种 方法 及 三 类 卫星 的 插值 精度 进行 研究 。 以 三 维 方式 较为 
形象 的 显示 了 两 种 插值 方式 下 ， 三 类 卫星 插值 精度 与 插值 阶 次 的 关系 ， 同 时 对 相同 阶 次 下 ， 
各 类 卫星 的 插值 精度 进行 了 统计 分 析 。 当 插值 阶 次 较 低 时 , 卫星 的 插值 精度 受 原始 精密 星 历 
的 影响 较 大 ， 会 呈现 出 同样 的 规律 性 ， 其 中 ，GEO 卫星 受到 的 影响 低 于 MEO 卫星 和 IGSO 
卫星 ， 随 着 插值 阶 次 的 增加 ,误差 的 规律 性 表现 会 逐渐 变 弱 。 三 类 卫星 采用 两 种 方法 插值 时 
的 最 佳 插值 阶 次 相差 均 在 1 以 内 ， 使 用 最 佳 插值 阶 次 插值 时 ， 滑 动 式 Lagrange 插值 精度 较 
非 滑 动 式 可 提高 11.96%~44.01%。 总 的 来 说 ， 三 类 卫星 中 ，GEO 卫星 在 相同 插值 阶 次 及 最 
佳 插值 阶 次 下 插值 精度 均 优 于 MEO 卫星 和 IGSO 卫星 ， 且 GEO 卫星 的 最 佳 插值 阶 次 也 较 
低 ; 两 种 方法 中 , 滑动 式 Lagrange 插值 精度 优 于 非 滑动 式 。 在 实际 使 用 时 ， 若 数据 量 较 多 ， 
滑动 式 Lagrange 插值 法 的 计算 量 会 明显 增 大 , 可 在 综合 考虑 插值 精度 和 插值 效率 的 情况 下 ， 
将 两 种 方法 组 合 使 用 。 
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